Threshold voltage and space charge in organic transistors 
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We investigate rubrene single-crystal field-effect transistors, whose stability and reproducibility 
are sufficient to measure systematically the shift in threshold voltage as a function of channel length 
and source-drain voltage. The shift is due to space-charge transferred from the contacts, and can be 
modeled quantitatively without free fitting parameters, using Poisson's equation, and by assuming 
that the density of states in rubrene is that of a conventional inorganic semiconductor. Our results 
demonstrate the consistency, at the quantitative level, of a variety of recent experiments on rubrene 
crystals, and show how the use of FET measurements can enable the determination of microscopic 
parameters (e.g., the effective mass of charge carriers). 



Organic single-crystal field-effect transistors (FETs) 
are opening new possibilities for the detailed investi- 
gation of the intrinsic electronic properties of organic 
semiconductors and of their interfaces [TJ |2]. Transis- 
tors where a single-crystal was suspended on top of a 
gate electrode, have led to the observation of intrinsic 
transport properties, such as mobility anisotropy [S] and 
metallic- like temperature dependence 4 . Current work 
is aiming at the systematic study of microscopic elec- 
tronic processes in these systems. Examples are the study 
of polaronic effects at the interface between organic crys- 
tals and highly polarizable dielectrics [5J [7] , the anal- 
ysis of band-like transport at interfaces with low-e ma- 
terials [H |9] , and the detailed investigation of electronic 
transport at metal/organic interfaces [10 . In most cases, 
a quantitative analysis of the data in terms of well-defined 
microscopic models has been possible, but the consis- 
tency of results obtained in different experiments remains 
to be verified. 

Virtually all experiments on single crystal FETs have 
focused on transport through a well-formed conducting 
channel, i.e. in the regime when the gate voltage is bi- 
ased well above the threshold voltage Vt- Here, we use 
rubrene single-crystal FETs for a systematic experimen- 
tal investigation of the behavior of the threshold volt- 
age itself. Specifically, we have measured the electrical 
characteristics of short channel transistors as a function 
of channel length L, and extracted the dependence of 
Vt on L and on source-drain bias Vds- We find that 
Vt systematically shifts to more positive values when L 
is decreased or Vds is increased, a behavior originating 
from changes in the space-charge transferred from the 
contacts into the semiconductor. We model the system 
using Poisson's equation, under the assumption that the 
density of states in the molecular crystals has the same 
functional dependence as in inorganic semiconductors, 
and find excellent quantitative agreement between exper- 
imental data and calculations, without introducing any 
adjustable parameter. Our results indicate that the effec- 
tive mass of carriers in the rubrene valence band is close 
to the free electron mass, provide information about the 



low-energy density of states, and show that the physical 
picture used to interpret a variety of recent experiments 
[PO"! ITT] is internally consistent at a quantitative level. 
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FIG. 1: Panels (a) and (b) show optical microscope images 
of devices used in this work. Panel (c) shows the transfer 
characteristics of FETs with different channel lengths, mea- 
sured at Vds — 60 mV (the curves are offset vertically for 
clarity). The sub-linear Va dependence is characteristic of 
contact dominated devices [10]. The vertical arrows point to 
the position of the threshold voltage in different devices. The 
two arrows with a continuous line indicate that for each data 
set the gate voltage is swept up and down. 

Not much is currently known about the microscopic 
processes determining the behavior of the threshold volt- 
age in organic transistors. Work done on thin- film de- 
vices has established that this behavior can be complex, 
exhibiting large sample-to-samplc fluctuations and insta- 
bilities [T2] . These phenomena have so far prevented sys- 
tematic experiments, and are posing severe technological 
problems. The situation appears to be different in single- 
crystal devices, with recent work showing considerably 
better reproducibility [13 . 

The fabrication of the rubrene FETs relies on the lam- 
ination of thin (< 1 /im-thick) single-crystals on top of a 
highly doped silicon substrate (acting as a gate) covered 
with 500 nm Si02, with predefined oxidized copper elec- 
trodes [TTj. The devices are identical to those that we 
used previously for the investigation of the contact resis- 
tance, and are therefore very well characterized [T0| 115] . 
In that work we focused on the high-gate voltage regime 
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with a fully formed channel; here we analyze the behavior 
of the threshold voltage Vt- 

Transport measurements were performed at room tem- 
perature, in the dark and in vacuum (2-7 1CT 7 mbar), 
using a HP4156A parameter analyzer. Owing to the 
short channel length, the resistance of the devices is 
dominated by the contacts, resulting in the character- 
istic (sub-linear) shape [10] of the transfer curves seen in 
Figs. 1 and 2. Special care was taken in restricting the 
applied gate voltage Vq to the narrowest possible range 
around Vt, to suppress possible effects of bias stress. Ac- 
cordingly, the output characteristics of our devices are 
hysteresis-free, with stable and reproducible Vt values. 
To extract these values, we have used the linear extrapo- 
lation method and fitted the linear regime of the transfer 
curve that follows the onset of conduction (as illustrated 
in the inset of Fig. 2). The uncertainty on the Vt values 
extracted using this method is at most ±0.1 V (for the 
shortest devices), and usually better. 

The behavior of the threshold voltage as a function of 
channel length L and bias Vds is apparent from Figs. 
[TJ; and [2] Fig. [TJ: shows the transfer characteristics of 
transistors with different channel lengths, measured at a 
same value of Vds- It is apparent that Vt systematically 
shifts to more positive values in devices with a shorter 
channel. Fig. [2] shows measurements done on an indi- 
vidual transistor for increasing values of Vds- A shift 
of Vt to positive values (see Fig. 3) is clearly present, 
which is found to be larger for shorter devices. These 
trends resemble those recently reported in Ref. [13] . Fig. 
[3] provides a complete overview of the behavior of Vt in 
our experiments. The different devices (i.e., sets of FETs 
with different channel length fabricated on a same crys- 
tal) investigated, exhibited identical trends and similar 
magnitude of the observed effects. 

The shift of Vt as a function of L and Vds origi- 
nates from the transfer of charge from the source and 
drain contacts into the semiconductor. At low Vds (i- e -> 
Vds — kT/e), charge transfer occurs to align the Fermi 
level in the metal and in the bulk of the semiconductor. 
As the contact separation is decreased, the correspond- 
ing charge density in the (bulk) region of the semicon- 
ductor between source and drain contact increases. At 
finite Vds, additional charge is transferred from the con- 
tacts to the semiconductor due to capacitive coupling, 
similarly to what happens in space-charge limited cur- 
rent measurements [16 . As charge density and potential 
are linked by Poisson's equation, the transferred charge 
induces a shift in electrostatic potential, which needs to 
be compensated by a shift in the threshold voltage to 
switch off conduction. In simpler terms, to switch off 
conduction, the space charge transferred from the elec- 
trodes into the bulk of the semiconductor needs to be 
compensated by an equal and opposite amount of charge 
accumulated in the channel (i.e., by a shift in Vt)- 

At a quantitative level we consider the case of low- 
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FIG. 2: Transfer characteristics of a rubrene single-crystal 
FET with channle length L — 600 nm, measured for different 
values of the source drain voltage Vds- The arrows indicate 
that in all measurements the gate voltage has been swept up 
and down. The inset illustrates how the threshold voltage is 
extracted from the data. 

bias first, and we calculate the profile of charge density 
transferred from the electrodes into the bulk of the semi- 
conductor (Fig. |4^i) by solving Poisson equation. To 
this end we employ a one-dimensional (ID) model, to de- 
scribe how electrostatic potential & s (x) and charge den- 
sity p s {x) vary in the bulk of the semiconductor, away 
from an interface with a metal contact (Fig. 4b). We 
then have 

<P*'( X ) = = g M H c-o*° M/kT . (1) 

dx 2 e e c eoec 

Here M H = 2(2nm p kTh~ 2 ) 3 / 2 , with m p the effective 
mass of holes in the "highest occupied molecular orbital" 
(HOMO) band of rubrene (i.e., the valence band of the 
organic semiconductor), eo is the vacuum permittivity 
and ec the relative dielectric constant of rubrene. The 
last equality assumes that the density of states in rubrene 
has the same dependence as in a conventional one-band 
semiconductor, i.e. parabolic dispersion relation with ef- 
fective mass m p , which allows expressing the local density 
of charge carriers p s {x) as a function of the local potential 
i> s (x) [T7]- Clearly, the validity of this hypothesis needs 
to be checked by comparing the results of the calculations 
to the experimental data. 

There are two boundary conditions associated to Pois- 
son's equation. The first, < & s (0) = $o/q, fixes the con- 
stant in the potential so that q$ s (x) corresponds to the 
local distance between the Fermi level and the valence 
band edge in the semiconductor (<S> is the energy dif- 
ference between the metal Fermi level and the valence 
band edge; when a large gate voltage is applied <&o cor- 
responds to the height of the Schottky barrier present 
at the metal/semiconductor interface). The second is 
d$ s /dx\ x = 0, which imposes that the electric field 



3 



6 



4 



> 



L = 3Wnm * 



4[3D nm 



500 nm 
600 nm 



2 




* 


■ 

* 






0- 


■ 

■ 

■ : t 


., i 

f * 


* 

■ i 


* 


600 nm 




• .»-« 
> # ■ 


■ * 


* * 


* * 


2 rii 

■ i 


-2 


* * * 


v - 


* ■ 




+ * 

3 |_JTl 




3 


0.5 


10 


1.5 


2.0 



FIG. 3: The dots represent Vr values of short-channel single- 
crystal devices as a function of Vi>s measured on devices with 
different channel length L (the dotted lines are guides for the 
eye). 



vanishes in the bulk of the semiconductor, far away from 
the interface with the metal. Under these conditions, the 
solution to Poisson's equation is 



2kT 



In 



' q 2 M H 

2kTe a e c 



(2) 



and the corresponding charge carrier density is given by 

-2 



p s (x) = qM H 



e *o/2kT + 



I q 2 M H 
2kTe e c 



(3) 



Note that in Eq. [I] we have neglected "bulk" contribu- 
tions to the charge density, e.g. due to (unintentionally 
present) dopants or thermally activated charge carriers. 
Accordingly, equations [2] and [3] are valid inside the semi- 
conductor, only at a small distance from the electrode 
surface, where the charge transferred from the contacts 
is larger than the one present due to dopants, which in 
high-purity crystals is N D m 2 10 14 crnT z [19 . This 
value is much smaller than what is usually found in or- 
ganic semiconductors, and hence in our devices the space 
charge penetrates deeper into the rubrene crystals. Us- 
ing Eq. [3] with $ = 0.13 - 0.15 eV [TO] (see discussion 
below) and taking into account the overlapping of the 
space charge regions coming from both contacts, it is 
easy to verify that the charge transferred from the con- 
tacts is larger than even for devices with L ~ 1 pm 
(normally, due to the large doping in organic semicon- 
ductors, this length is only a few tenths of nanometers 
[20]). 

To link the calculated quantities to the data, we impose 
that the total amount of space charge transferred from 
the contacts into the bulk region of the semiconductor 
between source and drain contacts is compensated, at 



threshold, by an equal and opposite amount of charge 
accumulated in the channel (i.e., by a shift in threshold 
voltage SVt(L)). We then get 



SVt(L) 



qt c < p s {L) > 
C G 

2qt c M H 



Le 



/ q 2 M H 
2kTe ec 



(4) 



where < p s (L) > is the spatial average of the carrier 
density present in the overlapping space-charge regions 
originating from the two contacts, Cq is the gate capaci- 
tance per unit area, and tc is the thickness of the rubrene 
crystal (measured with a profilometer) . This estimate of 
SVt is approximate, in that our calculation does not take 
into account the precise geometry. For our devices, we 
estimate that the error is a geometrical factor of the order 
of unity (see the below discussion on the Vds dependence 
of SV T ). 
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FIG. 4: Panel (a) schematics of the region between the source 
and drain contacts of a FET device. The arrows indicate 
charge transfer from the contacts into the semiconductor; 
the transferred charge is represented by the shaded region. 
Panel (b) shows the potential energy band-diagram near one 
metal/rubrene interface. Panel (c) shows the comparison be- 
tween the measured SVt{L) and the values calculated using 
Poisson's equation. The blue dotted line and the red line 
represent the values calculated with <E>o = 0.13 and 0.15 eV, 
respectively (corresponding to the known range of measured 
values for the CuO^ /rubrene interface). The inset shows that 
d[SVr(VDs)]/dVDS scales linearly with L~ 2 . 

Fig. HJ; shows the comparison between the predictions 
of equation [4] and the measured values 5Vt(L) = Vr(L) 
— < Vr > (where < Vr > = —1-9 V is the average 
threshold voltage obtained from devices with 10 pm < 
L < 50 pm; all measurements were performed at Vds — 
60 mV). The lines represent the results of the calcula- 
tions using the values of <& for a CuO^/rubrene interface, 
known from our study of the contact resistance: <&o typi- 
cally varies between 0.13 and 0.15 eV in different devices 
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10J. For the hole effective mass m p , we take the value 
extracted from recent optical spectroscopy experiments, 
which is close to the free electron mass All other 

quantities in the theoretical expression are known, and 
there are no free adjustable parameters. As it is appar- 
ent from Fig. [Ip, the quantitative agreement between 
calculations and experimental data is remarkably good. 

This result has several important implications. First, 
the comparison shown in Fig. [4j; relics on parameters 
extracted from completely different experiments, such as 
measurements of contact resistance and infrared spec- 
troscopy. Therefore, the agreement found in our analysis 
of SVrh(L) indicates that our description of the electronic 
properties of rubrene is internally consistent at a quanti- 
tative level. Second, our result indicates that describing 
the low energy density of states in the valence band of 
rubrene in terms of the "textbook" expression for con- 
ventional inorganic semiconductors is a good approxima- 
tion. Third, and more in general, this work shows how 
"simple" measurements of FET electrical characteristics 
can be used to extract microscopic parameters — e.g., the 
carrier effective mass — that are not otherwise easily ac- 
cessible experimentally. 

The analysis of the threshold voltage behavior as a 
function of source-drain bias confirms that the shift in Vr 
is due to space charge injected from the contact into the 
semiconductor. In this case we estimate the amount of 
charge injected into the semiconductor as CcVdSi where 
Cc is the capacitance between the injecting electrode and 
the rubrene crystal (this concept is identical to that used 
in the description of space-charge limited current I — V 
curves HI]). By reasoning analogously to the case of the 
L-dependence, we then obtain 



SV T (V DS ) 



Cc 
C G 



DS 



ec tctD 



DS 



(5) 



Where for the capacitance between contact and crystal, 
we have taken as a very simple approximation the ex- 
pression of a parallel plate capacitance. 

Equation [5] predicts that the shift of threshold volt- 
age is linear in source-drain voltage, which is the case 
for Vds > kT/q (see Fig. |3|, and that the slope of 
this linear relation should scale with the separation be- 
tween source and drain as £~ 2 , as found experimentally 
(see inset of Fig. |2J;). From the linear relation be- 
tween dl5Vr(VDs)]/dVDS and L~ 2 we extract the value 
of ectctD /e_D which is comparable to that obtained by 
directly measuring the crystal thickness (deviations of a 
factor of 2 — 3 are found in different samples, as can be 
expected given the crude estimates of the capacitances) . 

Note that all the effects that we have investigated here 
scale inversely to the gate capacitance Cq (see Eq. [4] 
and [5} . This is a typical signature of the so-called short- 
channel effects, which are well-known in silicon devices 
|21j . What is remarkable in our organic transistors is 
that, owing to the low doping concentration, these effects 



dominate already at a fairly large channel length (L ~ 
1 — 2 11m), i.e. length scales which are typical of devices 
used in practical applications [22] . 

In summary, rubrene single-crystal FETs are suffi- 
ciently stable and reproducible to perform systematic in- 
vestigations of threshold voltage shift. From the length 
dependence of this shift, we have extracted information 
about the microscopic properties of rubrene, such as den- 
sity of states and a quantitative estimate for the effective 
mass in the valence band. In conjunction with a variety 
of earlier experiments — infrared spectroscopy, and quan- 
titative studies of bias and temperature-dependent con- 
tact resistance — our results show that our current picture 
for the understanding of organic semiconductors does ac- 
count for many experimental observations in a way that 
is internally consistent at a quantitative level. 
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